Magneto-elastic effects and magnetization plateaus in two dimensional systems 



Samuele Bissola, 1 Valeria Lante, 1 Alberto Parola, 1 and Federico Becca 2 
1 Dipartimento di Fisica e Matematica, Universita dell'Insubna, Via Valleggio 11, 1-22100 Como, Italy 
2 CNR-INFM Democritos National Simulation Center and International 
School for Advanced Studies (SISSA), via Beirut 2-4, 1-34014 Trieste, Italy 
(Dated: February 4, 2008) 

We show the importance of both strong frustration and spin-lattice coupling for the stabiliza- 
tion of magnetization plateaus in translationally invariant two-dimensional systems. We consider 
a frustrated spin- 1/2 Heisenberg model coupled to adiabatic phonons under an external magnetic 
field. At zero magnetization, simple structures with two or at most four spins per unit cell are 
stabilized, forming dimers or 2 x 2 plaquettes, respectively. A much richer scenario is found in the 
case of magnetization m = 1/2, where larger unit cells are formed with non-trivial spin textures 
and an analogy with the corresponding classical Ising model is detectable. Specific predictions on 
lattice distortions and local spin values can be directly measured by X-rays and Nuclear Magnetic 
Resonance experiments. 
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A magnetization plateau occurs when the magnetiza- 
tion remains constant over a range Ah of applied mag- 
netic fields. The width of the plateau can be expressed 
in term of the excitation spectrum Ah = E(M + 1) — 
2E(M) + E{M - 1), where E{M) is the total energy at 
fixed magnetization M (measured in units of g/is). This, 
in turn, implies that the energy (per site) as a function of 
the magnetization (per site) displays a cusp-like singular- 
ity in the thermodynamic limit. In general, plateaus are 
absent in classical models with magnetic ground states 
whenever the magnetization is not collinear with the 
field. [H Instead, they identify particularly stable quan- 
tum phases characterized by a spin gap. 

Usually, gaps in the excitation spectrum directly reflect 
the structure of the primitive cell of the lattice according 
to the commensurability condition [3] 
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where S is the magnitude of the spin, m the magneti- 
zation per site in units of gfisS (i.e., the average on- 
site spin component parallel to the magnetic field), and 
£ the number of spins in the primitive cell. According to 
Eq. {I}, magnetization plateaus in models with a single 
spin-1/2 per unit cell may occur only as a consequence 
of the spontaneous breaking of the translational symme- 
try, leading to an £-io\d degenerate ground state. Al- 
though this relation has been rigorously proved in one- 
dimensional (ID) systems 0|, it is believed to have a 
much wider validity. [3] In fact, several analytical and 
numerical calculations in chains and ladders have defini- 
tively confirmed the presence of various magnetization 
plateaus at the expected positions. 0, [H, @ Moreover, 
few important cases of magnetization plateaus in higher 
dimensions are also known. The most relevant exam- 
ple is given by the orthogonal dimer Heisenberg model, 
which closely represents the structure of copper planes 
m SrCu 2 (B0 3 ) 2 . [7] In this case, the theoretical predic- 
tions have been confirmed by experimental evidence of 



magnetization plateaus at m = 0, 1/8, and 1/4. [8j Some 
evidence for a m — 1/3 plateau has been proposed for the 
triangular lattice, [l(| while in the square lattice, the 
m = properties of the J\ — Ji model are still debated. 
Indeed, although for J2/J1 ~ 1/2 the ground state is be- 
lieved to be disordered, the existence of a finite triplet 
gap is much less clear, [Til, E3] casting some doubt as to 
the possibility of having an m — plateau. 

The interest in the J\ — Ji model has grown due to 
the recent discovery of two materials well described by 
a two-dimensional (2D) quantum antiferromagnet, i.e., 
Li 2 VOSi0 4 and VOM0O4. [HEi Although the experi- 
mental magnetization curve of these compounds has not 
yet been considered, the magnetization properties of the 
J\ — Ji model have been recently considered by using ex- 
act diagonalization calculations, leading to some evidence 
in favor of a magnetization plateau at m = 1/2. [l5| This 
outcome has been interpreted as a consequence of the 
emergence of a 2 x 2 super-cell. According to this sce- 
nario, inside each cell the magnetic moments acquire a 
preferential orientation along the direction of the mag- 
netic field, leading to a configuration with three up and 
one down spin. [151 ] Unfortunately, the presence of such 
a half-magnetization plateau is limited to a very nar- 
row region close to J2/J1 ~ 1/2, revealing the difficulty 
of stabilizing such a state in the pure spin model with- 
out other degrees of freedom. In this respect, the spin- 
lattice coupling represents one of the most relevant phys- 
ical mechanisms to enhance the stability of phases break- 
ing some of the lattice symmetries, making it easier to 
identify these states in numerical studies of small lat- 
tices. From general arguments, the super-exchange cou- 
plings are ultimately generated by the virtual hopping of 
electrons through neighboring sites and strongly depend 
upon lattice distortions. [l6| The role of the spin-lattice 
coupling in frustrated spin systems has been extensively 
considered in the absence of an external magnetic field 
for both ID [13, OH and 2D systems. [H| More recently, 
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the importance of lattice distortions for stabilizing mag- 
netization plateaus has been discussed in a simple ID 
spin model, [2(| and in a classical Heisenberg model on 
a pyrochlore lattice. 21 1 

In this Letter we investigate the possible occurrence of 
magnetization plateaus in the frustrated spin- 1/2 Heisen- 
berg model on the square lattice by the numerical anal- 
ysis of the periodicity of the distortion pattern induced 
by spin-lattice coupling. The model is defined by 
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where is the spin-1/2 operator at site i, 5vi is the dis- 
placement of atom i, assumed to be in the plane, and 
d%j = ||r< — TjH is the distance between atoms i and j. 
The sum runs over nearest [i.e., J\ = .7(1)], second [i.e., 
Ji = J(\/2~)] and third [i.e., J3 = J(2)] neighbor sites on 
a 2D square lattice, while c?°- = ||R° — R§|| is the distance 
between sites i and j in the undistorted lattice. Energy 
units are fixed by the choice J± = 1 . In order to limit the 
number of parameters, only the elastic constant between 
nearest-neighbor sites [K = K(l)] has been included. 
However, we checked that our conclusions are not qual- 
itatively affected by considering further neighbor elastic 
constants. In transition metal compounds, the super- 
exchange theory combined with empirical dependences 
of hopping integrals on distance [16[ leads to exchange 
integrals that vary like the inverse of the distance to a 
given power a M (/x = 1,2,3). For small displacements we 
can write: 
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with Sdij = dij - d% ~ (R° - H°) • ($r< - 5r 3 )/d%. 
Since the Hamiltonian is invariant under the rescaling 
a M = Xa^, K = X 2 K and 5rj = Sri/X, K can be fixed. To 
study this system, we adopt the Lanczos diagonalization 
technique in finite clusters with periodic boundary condi- 
tions, which allows for an unbiased determination of the 
lowest-energy configuration at fixed distortion pattern. 
In order to determine the optimal configuration of lat- 
tice displacements <5r, we use an iterative procedure. [IH 
Since the local displacement of each spin from its equi- 
librium position is considered, we are able to describe 
all kinds of distortions, including the dilation and the 
shrinking of the lattice. We performed systematic calcu- 
lations at different frustrations and lattice couplings in a 
4x4 cluster at m — and m = 1/2, i.e., at the magneti- 
zation values where plateaus are commonly expected in 
this geometry. 

At vanishing magnetization, the classical J1—J2—J3 
model on a rigid lattice shows four magnetically ordered 
states: Neel with momentum (7T,7r), collinear with mo- 
menta (tt, 0) and (0, tt) and two helicoidal phases with 
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FIG. 1: Ground-state phase diagram of the Ji — J2 — Jz model 
for m = 0, K = 10, and ,h = (a) or ,h = 0.5 (b). See text 
for the precise description of the various phases. 



momenta (q, q) and (q, 7r), (tt, g), where q varies continu- 
ously with the parameters in the Hamiltonian. [22, 2^, 24 1 
Quantum fluctuations, enhanced by the competing an- 
tiferromagnetic couplings, can drive the system away 
from these semiclassical behaviors and stabilize uncon- 
ventional quantum phases without long-range magnetic 
order, which may eventually display a magnetization 
plateau. At J3 = 0, the magnetically disordered regime 
is widely believed to occur around the maximally frus- 
trated point (i.e., J2 = 1/2), while at J2 = a non- 
classical phase appears between the Neel and the spiral 
state (q, q) close to J3 ~ 1/2. However, in both cases, the 
nature of the disordered phases is still controversial and 
several proposals appear in the literature: 25j Valence 
bond crystal (VBC) columnar states 0, [i3] or spin liq- 
uids. 11, Ii| Very recently a singlet state with plaquette 
order has been claimed to be stabilized around the line 
J 2 + J 3 ~ Ji/2 with J 3 > 0. [2^ 

When the spin-lattice coupling is taken into account, 
we expect that each of the previous spatial orderings 
would lead to a characteristic distortion pattern which 
can be easily identified. In Fig. [1] we present the results 
for the ground-state phase diagram obtained by Lanczos 
diagonalizations of the Hamiltonian J2} in the J^-a plane 
for = a and K — 10. Two values of J2 — and 1/2, 
representative of the weak and strong frustration regime, 
will be considered. 

Among all possible lattice deformations, only four are 
stable at m = 0, depending on the values of the parame- 
ters (see Fig.[T]): i) The square lattice with uniform bond 
lengths: All spatial symmetries are preserved and the 
presence of the spin-phonon coupling a just leads to a 
renormalization of the bond lengths, ii) The dimerized 
lattice with two different bonds in one direction and one 
bond in the other one. This phase breaks both ir/2 rota- 
tion and the translational symmetry along one direction 
and the ground state is four-fold degenerate. In this case, 
the commensurability condition {1} is satisfied leading to 
spin gap and a magnetization plateau. Hi) The square- 



3 



plaquette phase, with dimerization in both directions. In 
this state, rotational symmetry is preserved since a 2 x 2 
primitive cell is stabilized. Also in this case the ground 
state is four-fold degenerate and the spin gap is finite. 
In any case, the plausible occurrence of spiral magnetic 
ordering with incommensurate q at large J3 suggests to 
be cautious with the numerical results in small lattices 
f° r J3 ^ 0.6. For J 2 = 1/2, we also found iv) the rect- 
angular phase with different bond lengths in the x and 
y directions. This phase breaks the 7r/2 rotational sym- 
metry but is translationally invariant, as expected when 
collinear magnetic ordering is present. In this case the 
ground state is two-fold degenerate, leading to a finite- 
temperature phase transition. [30I ] Finally, the absence 
of distortions for low values of J3 is consistent with a 
magnetic ground state displaying Neel order. In all of 
these cases, the local magnetization (Sf) is uniform in 
the whole lattice. 

In order to better understand the physical properties 
of the two intermediate phases, ii) and Hi), we have an- 
alyzed the energy spectrum both in the 4x4 and in the 
6x6 rigid lattice. Some evidence of translational sym- 
metry breaking emerges from the quasi degeneracy of the 
states with momenta (0,0), (0, 7r) and (tt,0). However, 
the ordering of energy levels on these small systems is 
not able to discriminate between a columnar dimer or a 
plaquette state, because both the d-wave, zero momen- 
tum state and the (71", tt) s-wave singlet have comparable 
energies in these clusters. 

Let us now move to the more interesting case of m = 
1/2. The phase diagram is shown in Fig.[2j At small J3, 
a semiclassical scenario is consistent with both the shape 
of the lattice and the observed uniform local magneti- 
zation (i.e., (Sf) — 1/4). Indeed, an applied magnetic 
field cants the spins which preserve a magnetic ordering 
in the plane orthogonal to the field direction. The nu- 
merical evaluation of the spin-spin correlation function 
by means of Lanczos diagonalizations on the 4x4 and 
6x6 rigid cluster is fully compatible with the semiclas- 
sical scenario both at small and large J2, namely weak 
correlations along the magnetic field and marked peaks of 
the structure factor in the orthogonal plane. The analyt- 
ical evaluation of the susceptibility via spin-wave theory 
shows that, for weak phonon coupling, the lattice un- 
dergoes global deformations but bond lengths and local 
magnetizations remain uniform throughout the lattice, 
thereby inhibiting magnetization plateaus. 

By increasing the third neighbor coupling J3, in ad- 
dition to the same lattice deformations already found at 
zero magnetization, other two phases emerge (see FigJ5]): 
v) The "trapezoidal" phase characterized by a 4 x 2 primi- 
tive cell made of two aligned congruent isosceles trapezia, 
tilted by tt one respect to the other. In this case, the lat- 
tice dilates in one direction and exhibits three different 
bond lengths, still preserving a rectangular shape (see 
Fig. E]) . The local magnetization is not uniform but it ac- 
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FIG. 2: The same as in Fig. [T] but for m = 1/2. 
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FIG. 3: Distortion pattern of the Ji— J2 — J3 Heisenberg 
model on the 4x4 lattice at m = 1/2 in the "trapezoidal" 
phase at J2 = 0.5, J3 = 0.2, and a = 3.5 (a) and in the "clas- 
sical Ising" phase at J2 = 0.5, J3 = 0.4, and a — vT(5 (b). 
Distortion pattern of the Ising model on the 8x8 cluster at 
J2 = 0.5, J3 = 0.4, and a — \/l0 (c). Stars mean (Sf) ~ 0, 
while empty and full dots indicate positive and negative lo- 
cal magnetization close or equal to the maximal value 1/2, 
respectively. 



quires two different values (see Fig. [4j . vi) The "classical 
Ising" phase characterized by a 4 x 4 unitary cell com- 
posed of four identical scalene trapezia which can be ob- 
tained one from the other by means of suitable rotations 
(see Fig. [3]). The lattice is not squeezed in any direction, 
but develops a complex bond pattern with eight different 
bond lengths and three values of the local magnetization: 
Two positive and one negative, with modulus close to 
the limiting value S — 1/2 (see Fig. 0]). This rather rich 
texture of local magnetizations (Sf) is stabilized for non- 
zero third neighbors coupling and can be directly tested 
by Nuclear Magnetic Resonance (NMR) experiments. In- 
terestingly, this spin pattern has a direct interpretation 
on the basis of a purely classical Ising model. Indeed, we 
checked that the complex topology previously described 
is exactly shared by a 2D Ising model with antiferromag- 
netic interactions up to third neighbors and spin-lattice 
coupling at one-half of the saturation magnetization. A 
numerical study of this model in the same parameter 
regime has been carried out both in the 4x4 and in 
the 8x8 lattices with similar results: The spin patterns 
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Li 2 VOSi0 4 or VOM0O4. 

We acknowledge interesting discussion with P. Car- 
retta and D. Poilblanc. This work has been partially 
supported by COFIN 2005 and CNR-INFM. 



FIG. 4: Behavior of local magnetization (SI) as a function of 
J3 at m = 1/2 and a = \/T0 for J2 = (upper panel) and 
J2 = 0.5 (lower panel). 



is characterized by the requirement of placing the down 
spins so to avoid first and third neighbors, while allow- 
ing a single second neighbor down spin, see Fig. [U [3l| 
It should be noticed that this constraint leads to slightly 
different patterns in the 4x4 and the 8x8 clusters, be- 
cause of periodic boundary conditions, and on the 4x4 
is compatible with the results on the quantum Heisen- 
berg model. Moreover, also the distortion pattern of the 
classical model is extremely similar to that obtained in 
the quantum model (see Fig.[3|). 

It is worth noting that once the magnetic energy gain 
overcomes the elastic energy, the transition lines among 
the different phases depend only weakly on the spin- 
phonon coupling a. This suggests that, in the thermo- 
dynamic limit, the distorted phases might remain stable 
up to arbitrarily small a, as can be analytically shown 
for the nearest-neighbor Ising model. In such a case, the 
m = 1/2 magnetization plateau would be present, in the 
J1 — J2 — J3 model, independently of the strength of spin- 
lattice coupling. 

In summary, we showed that magnetization plateaus 
can be stabilized in the square lattice Ji~J2~J3 anti- 
ferromagnetic Heisenberg model in the presence of spin- 
phonon coupling. Third neighbor coupling appears to 
be an essential ingredient for the appearance of plateaus 
both at zero and one-half magnetization. While at m = 
the lattice dimerizes as expected on the basis of a VBC 
phase, at m — 1/2 a novel and complex distortion pat- 
tern characterizes the broken symmetry state. This phase 
can be faithfully interpreted in terms of a frustrated 
Ising model coupled to the lattice and appears as a gen- 
uine result not much affected by the finite size of the 
cluster we analyzed. The emerging scenario is consid- 
erably richer than predicted in a seminal study of the 
J\ — J2 model on a rigid lattice [l5| and identifies a sim- 
ple mechanism for the stabilization of complex spatial 
structures. The on-site magnetization displays a distinc- 
tive pattern which may be measured by NMR experi- 
ments in strongly frustrated 2D magnetic materials, like 
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